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Abstract 

Weak gravitational, electromagnetic, neutrino and scalar fields, considered as 
perturbations on Kerr background satisfy Teukolsky Master Equation. The two 
non-trivial equations obtained after separating the variables are the polar angle 
equation and the radial equation. We solve them by transforming each one into 
the form of a confluent Heun equation. The transformation depends on a set 
of parameters, which can be chosen in a such a way, so the resulting angular 
and radial equations separately have simple polynomial solutions for neutrino, 
electromagnetic, and gravitational perturbations, provided some additional con- 
ditions are satisfied. Remarkably there exists a class of solutions for which these 
additional conditions are the same for both the angular and the radial equations 
for spins |s| = 1/2 and |s| = 1. As a result the additional conditions fix the de- 
pendence of the separation constant on the angular frequency but the frequency 
itself remains unconstrained and belongs to a continuous spectrum. 

Keywords: Teukolsky master equation, Heun equation, Heun functions, perturbations 
of black holes. 



1 Introduction 



Fields of various types - scalar, neutrino, electromagnetic, and gravitational - have 
been extensively studied as perturbations to known solutions of Einstein's equations 
for configurations with spherical and cylindrical symmetry. The fields considered are 
weak in the sense that we can neglect the influence of their stress-energy tensor on the 
background metric. Regge and Wheeler jl] and Zerilli [21 13] were the first to study the 
linear response of Schwarzschild solution of Einstein's equations to perturbations. In 
order to study Kerr metric perturbations Teukolsky |1]-[T] analyzed the components 
of Weyl tensor, using Newman- Penrose formalism [11]. (For a detailed extended pre- 
sentation see [I2].). As a result one obtains the Teukolsky Master Equation, which 
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describes the dynamics of various fields of different spins as perturbations to Kerr 
metric. In recent years there is an increased interest on the subject [T3]-[16j, mostly 
aimed at studying the quasi-normal modes. Another problem analyzed via Teukolsky's 
equations is related to the decaying of the various fields present during a gravitational 
collapse at very late times at large distances - the so called late-time tails. All these 
investigations however, are performed via indirect, approximate methods [T3t [T^.[T6]- 

m- 

On the other hand, already for quite some time it has been recognized in the litera- 
ture [ini Eni EH 122] that Regge- Wheeler and Teukolsky's equations can be transformed 
into the form of a confluent Heun equation [23]- [29]. The reason there has not been 
much attention paid to the Heun-type solutions is that they are not completely ana- 
lyzed and, in general, difficult to work with. Some basic classes of exact solutions to 
Rege- Wheeler equation in terms of special solutions to the confluent Heun equation - 
the so-called confluent Heun functions (see the Appendix), were described recently and 
were used for finding solutions to a number of physical problems [5D]-[5^. 

Following the articles [30]- [32] we continue with the application of the confluent 
Heun functions to Teukolsky's equations. The first results, presented in [33]- [38], were 
very encouraging and drew special attention to the solutions in terms of the confluent 
Heun polynomials [23]- [29] (see the Appendix). It should be emphasized that long time 
ago in [ini [20] it was recognized by Baldin, Pons, and Marcilhacy that the conditions 
for polynomial solutions to Heun equations lead to polynomial solutions to Teukolsky's 
equations in a generalized sense, i.e. polynomials multiplied by non-polynomial factors 
which are elementary functions. Having in mind the general description of all 256 
classes of factorized solutions to Teukolsky Master Equation [391 HQ] we intend to 
focus on the mathematical properties of some of them and study various physical 
applications. 

The general description of all polynomial solutions of Teukolsky Master Equation 
was given for the first time in [39| HQ] . These fall into two different classes. For the 
first class, the first polynomial condition (110. 9p . called the ^Ar-condition in [29], [391 HQ], 
is automatically satisfied. For waves of spin |s| this condition fixes only the degree 
(A^ + 1) = 2|s| of the second polynomial condition Atv+i = 0. For the second class of 
polynomial solutions the iJat- condition is fulfilled only for certain complex frequencies 
un which belong to definite equidistant discrete spectra. For the two classes the second 
polynomial condition Ajv+i = defines an algebraic equation of degree 2|s| for the 
second separation constant: Em = E,m{<^), and Em = Em{ujN), correspondingly. 

Here we are considering only polynomial solutions of the first class. Thus an inde- 
pendent derivation of the specific relations, valid only for the first class of polynomial 
solutions becomes possible. It is based on a direct check of the two necessary conditions 
(110.91) and (110.101) . which together are sufficient to ensure the polynomial character of 
the solutions (See the Appendix.). 

Below we present an independent derivation of the first class polynomial solution 
both for Teukolsky's angular and radial equations using the notations of reference [28] . 
This notation has some advantages since it simplifies significantly the form of the 5n- 
condition. The correspondence between the notation of [28] and the notation used in 
[251126], [29],[3Q]-[lQ] and in the computer apphcation Maple is described in section 9.4 



2 



of the Appendix. 

In the present paper we consider a specific type of evolution of week fields with 
spin |s| = 1/2,1 and 2 on Kerr background. The solutions studied here are double 
polynomial solutions that describe one-way waves of corresponding spins, the so-called 
total transmission modes. These are factorized solutions to Teukolsky Master Equation, 
in which the solutions both of the angular and the radial equations (of the same spin 
weight s) belong to the corresponding first classes of polynomial solutions, introduced 
in [39lB0] . Here we show that these solutions yield a complex one-parameter continuous 
spectrum of the frequency u and derive the explicit form of the separation constant 
E in the various cases. Finally we discus some overall solutions of Teukolsky Master 
Equation, constructed making use only of these continuous spectrum solutions. 

To the best of our knowledge this is the first time when for Teukolsky Master 
Equation exact solutions with continuous spectrum are presented for a specific bound- 
ary problem. An interesting observation is that continuous spectrum emerges only 
for neutrino and electromagnetic waves, because of the simultaneous fulfillment of the 
polynomial condition both for the angular and the radial Teukolsky equations. We 
have to stress that such a simultaneous fulfillment is not in place for gravitational 
waves. The physical consequences of this mathematical result may be deep and very 
important. Its roots can be traced back to some results, originally obtained in [7] and 
developed further in [12]. We present here the mathematical basis, needed for further 
developments in this direction. 

In the next section we start by reminding the procedure for separation of the vari- 
ables in Teukolsky Master Equation via factorization of the solutions and the corre- 
sponding basic results. In section 3 we present the general scheme for transforming 
Teukolsky's radial equation (TRE) into the one of the many known "canonical" forms 
of the confluent Heun equation p5]- |29j^ . namely into the so-called non-symmetrical 
canonical form. We show that for specific values of the indices of the regular singular 
points [41j the first condition for polynomial solution to TRE in the form of a con- 
fluent Heun equation is automatically attained. We impose the second condition for 
having a polynomial solution and obtain the value of the separation parameter E as 
function of the frequency u. In section 4 we continue by presenting the transformation 
to non-symmetrical canonical Heun form of Teukolky's angular equation (TAE) and 
again show that for corresponding specific choice of the indices of its regular singular 
points we obtain again a polynomial solution. In order to achieve this result we derive 
the explicit form of the second polynomial condition and arrive at the result that in 
some cases it is the same as for the radial equation. Thus we find a simultaneous 
fulfillment of the polynomial conditions for the angular and for the radial Teukolsky's 
equations for perturbations with spin \s\ = 1/2 and 1. In section 5 we discuss why such 
a simultaneous fulfillment of the polynomial conditions is not possible for gravitational 
waves = 2). In section 6 we present the basic properties of the overall solutions 
to the Teukolsky Master Equation constructed only from the factorized solutions with 

^In the literature there still does not exist a commonly accepted standard form both of the five 
classes of Heun equations and of the corresponding Heun functions. Different possible forms are in 
use, since in different applications the authors prefer the form, which is most suitable for their specific 
needs. 
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continuous spectrum. We show that these solutions describe one-way colhmated waves, 
which may be regular along the rotational axes, despite the singular character of the 
polynomial solutions of the angular Teukolsky equation. In mathematical sense these 
solutions form a natural orthogonal basis of singular functions for integral representa- 
tion of physically meaningful solutions. In the conclusion we give a brief summary and 
ideas for future studies on the matter. 

In an Appendix some basic properties of the confluent Heun equation and its solu- 
tions and different forms are presented for the reader's convenience. 



2 Spin weight s fields on Kerr background 



In this section we present some basic results of Teukolsky's approach |l]-[7] to the 
perturbations of spin |s| of Kerr vacuum solution for the metric of a rotating black 
hole. In Boyer-Lindquist coordinates the metric is given by [12] , 



2Mr\ , , AaMrsin^e , , , E , , 

1 dt^ H dtdd) dr^ 

S / S ^ A 



- Ede^ - 



+ + 



2MaV sin^ 6 



sm 



(2.1) 



Here M is the Keplerian mass of the rotating black hole and a is its angular momentum 
per unit mass. Also, A and S are defined in the usual way: 



A = r - 2Mr + 



S = + cos^ 6 . 



(2.2) 



The dynamics of a massless field \E' = '^{t,r,6,(j)) with spin weight s is described by 
Teukolsky Master Equation: 
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(2.3) 



In the above equation we have the following expressions for 



For s = 1/2, \E' = Xo and for s = —1/2, \I' = p ^Xi, where Xo and Xi represent 
the two components of the neutrino spinor in Newman-Penrose formalism. 

For s = 1, \E' = and for s = —1, \I' = p~^(p2, where (po and are Maxwell 
tensor tetrad components in Newman-Penrose formalism. 

In the gravitational case for s = 2, = iJjq and for s = —2, = p~'^ilj4, where ipQ 
and are Weyl tensor Newman-Penrose tetrad components. 
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In all cases p = — l/(r — ia cos 6) and the choice of the Kinnersley tetrad is assumed. 
In order to separate the variables following [l]-[7] we set 



<f{t, r,e,(f)) = — Yl ^-(^' ^' ^' ^' ' (2-4) 

where \l/m(ci;; t, r, 6', 0) = e^'^'^e^^'^^ Sm{(^', 0)Rm{(^', r) , and e is a parameter defining the 
contour of integration. The azimuthal number has values m = 0, ±1, ±2, . . . for integer 
spin, or m = ±1/2, ±3/2, ... for half-integer spin HI [TOl HQ] . 

Thus we are looking for an integral representation with a factorizedkernel t, r, 6, 

For the unknown factors ^^.(w; 9) and Rm{^', r) we obtain respectively Teukolsky's an- 
gular equation (TAE): 

^ ^ '^sine^] S^{u;9)+ (2.5) 



sin V 96^ 

/ 2 2 2 /I n m'^ + s"^ + 2ms COS 9 ^ \ ^ , 

+ { a CO cos — 2acos cos ^ h Em Sm[u:;0) = , 

\ sin 9 J 

and Teukolsky's radial equation (TRE) 

A-^^ (a^+i|-) RUu^; r) ± ~ ~ ^ ± Aisujr - \^ Rm{uj- r) = , 

(2.6) 

where K = {r'^ -\- a?) uj — am and = Em + a^cu^ — 2amu — s{s + 1) is the separation 
"constant"[§. These are the two equations, some special solutions of which we will study 
in detail in our paper. 



3 Solutions to TRE 



3.1 Transforming TRE into the non-symmetrical canonical 
form of Heun equation 

Both equations (12. 5p and (12. 6 p are second order ordinary differential equations with 
two regular singular points at finite value of the independent variable and one irregular 
singularity at infinity. The most general equation with such properties is the confluent 
Heun equation [28]. Thus both TRE and TAE are specific cases of confluent Heun 
equations. At this point we continue by transforming the radial equation (12. 6p into 
one of the forms of the confluent Heun equation, the non-symmetrical canonical form 
(110. ip . In order to do so we apply the so-called s-homotopic transformation to (12.60 
by setting 

Rmiuj; r) = (r - r+)«(r - r_) Vi/(r) , (3.1) 

^In the literature often another form of the separatfon constant is used, namely A„i = £'.,„ — s(s + l) 
so A„i can also be written as A„i = Am + a^w^ — 2amuj. 

^This anzatz was used for the first time for analytical and numerical studies of the problem at 
hand by Leaver [18] , [43l |44] . 
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where r_|_ and r_ are the event and Cauchy horizons of a Kerr black hole defined 
hy r± = M ± \/ — and ^, ri (the indices of the regular singularities) and C, 
are parameters to be determined. By substituting (13. ip into (12.61) and after some 
straightforward algebra we arrive at the following equation for H{r): 

(pH{r) ^ + s + l 2r] + s + l \ dH{r) 

H I" 2Q 1 - h 



dr^ V r — r_i_ _ ^_ J dr 

1 



-{ [Au^M + 2C(^ + r] + s + l) + 2isuj] r + 



(r — r+)(r — r_) 

+ (e + V? + i^ + v) + 2s(e + r]) - 2C(er- + T/r+) - 2<M - 2CM - 

-2ac<;m + 4c<;2m2 - 2isujM - X„,}H{r) = . 

In order to obtain this form of the equation we had to fix the parameters ^, rj and (. For 
them we end up with quadratic equations and thus with pairs of possible expressions, 
namely: 

= ia+{uj — mQ+), (,2 = —s — ia+{uj — rriQ^) , (3.2) 



rji = —s + ia_{u! — mfl^), ri2 = —ia_{u — mfl^) , (3.3) 



and 



Ci = i^, C2 = , (3.4) 

where we have set 

2Mr± a 
r_|_ — r_ 2Mr-|- 

Any one of the eight possible triplets of expressions from (13. 2p . (13.31) . and (13.41) leads to 
a Heun equation in the desired form with different parameters. In order to determine 
these parameters it is necessary to perform one more step and to introduce dimension- 
less variables instead of r. In order to keep the symmetries in the problem manifest it 
is best to set different variables for positive and for negative spin weights, namely for 
s = 1/2, 1, 2 and for s = —1/2,-1, —2 we will have respectively: 

r_(_ — r J- — j'_ 
^= , -z= . 



r+ — r_ 



Note that from now on in the paper we will denote the sign of the spin weight (and 
when there might be a confusion the spin weight itself) by a subscript to the left of the 
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variable. Also always the signs "+", or ± used as subscripts to the right from the 
variables will be related to the two singularities in the radial equation, namely to the 
event and the Cauchy horizons. After performing these adjustments we end up with a 
Heun equation in the non-symmetrical standard form ( 110. 1^ : 

— — + 4(±p) + — + — — - + — --- — H = (3.6) 

with the parameters given by the following expressions: ±p = =f(^+ ~ ''"-)C/2, +7 = 
_5 = 2^ + s + 1, _7 = +5 = 27] + s + 1, 



±a = 2Moj^C^ + ^ + r] + s + l + isujC-\ (3.7) 



and 



±a = -2Cr± + ^ + V + s + l + —j - [(^ + + + + + v) - 



-2C({r_ + r/r+) - 2sCM - 2CM - 2aujm + 4w^M^ - 2isujM - \m] , 

for any choice of a triplet ^, 1], (. In general there exist pairs of solutions to Heun 
equation expressed as appropriate power series about each one of the singularities. The 
solutions of the equation (13. 6p can be expressed with the use of the so-called Frobenius 
solutions Hc^^\±p, ±a, ±7, ±S, ±a; ±z). There exists a multitude of other solutions to 
Heun equation which can be obtained from these Frobenius solution by interchanging 
the finite singular points or performing appropriate s-homotopic transformations (See 
the Apendix). 



3.2 Polynomial solutions 

Thus we have completed the transformation of TRE to the non-symmetrical form of the 
confluent Heun equation and we can formally identify the solutions which are given by 
specific Heun functions. There exists though a special case, in which the confluent Heun 
equation admits polynomial solutions (See the Appendix). This special case depends 
on the values of ±a in (13. 7p and another, more involved condition on the parameters 
in the Heun equation. If the parameters ±a are equal to negative integer numbers or 
zero and if we can solve the second condition then we will have polynomial solutions 
to Heun equation. At this point we can make the following observation regarding the 
possible values of ±a in (13. 7p : In the cases of positive spin weights if we pick the 
values ^2, Viy (2 from ( 13. 2p . (13. 3p . and (13. 4p . then we obtain = 1 — 2s, which 
equals for neutrino, —1 in the electromagnetic case and —3 in the gravitational case. 
Similarly, for s = —1/2, s = —1, or s = —2 we have to choose ,^1, 772, and Ci and 
we will get _a = 1 + 2s, which again gives a = 0, — 1 or —3 for the neutrino, for 
the electromagnetic, and for gravitational case respectively. Thus we have identified 
particular combinations of the values of the parameters ^, r], and ( for which we may 
expect to find polynomial solutions of the radial equation for |s| = 1/2, 1,2 since the 
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necessary condition a = — A^, A^-integer, is satisfied. Also, it can be seen easily that for 
s = we obtain = 1 so there is no polynomial solution for scalar fields. The second 
condition will depend on the value of A^. It is an algebraic equation of order A^ + 1 
and leads to polynomial solutions to Heun equations of order A^. Thus we expect that 
in the neutrino case the polynomial solutions, if they exist, are simply constants, in 
the electromagnetic case - linear functions, and in the gravitational case the solutions 
are cubic polynomials. We will consider here in more detail only the electromagnetic 
case. The results for s = ±1/2 can be easily obtained using the same procedure. The 
polynomial solutions to TRE with spin |s| = 2 were described long time ago in quite a 
different setting by Chandrasekhar |45j. For a more recent treatment see ^6] . 



3.3 Polynomial solution to TRE for spin \s\ = 1 

In order to find a polynomial solution to TRE we have to impose in addition to the fact 
that ±a = —N the condition fllO.lOp . In this case we are looking for an expansion about 
the singular point at infinity. The second polynomial condition for electromagnetic 
perturbations has the form g^^^^g^^ = h^i^ fl[^ with the coefficients from the three-term 
relation (110.81) given by ^^q^^ = —a — 4p + 7 + 5, g^^ = —a, flf^ = —4:p, hf^ = —7. 
Provided we make the above mentioned choices rendering ±a = —1 we obtain the 
following values for the remaining parameters in the Heun equation: ±p = {i/2)uj{r+ — 
r_),±7 = ±2za_(a; — mQ-),±6 = ^2ia^{uj — mQ+),±a = ±Em + a^u;^ — 2aujm =F 
2iu:r±. This second polynomiality condition essentially fixes ±a. It leads to a quadratic 
equation for ±a with solutions ±0"! = ^2iu}r± + 2^J au}{au — m) and ±a2 = T'^ioor^ — 
2-\J auj{abj — m) (recall that the ± sign to the left of a refers to the sign of the spin 
weight, while the subscripts 1 and 2 number the solutions to the quadratic equation). 

The first important result from solving the second condition for having polynomial 
solutions is that we find the dependence of the separation "constant" Em on the fre- 
quency u, which is the same for both s = +1 and for s = — 1 [33 SO] and is given 
by 

±iEm{auj)i = —a^iJ^ + 2aum + 2 a/ au^au — m) 
±iEjn{auj)2 = —a^uj'^ + 2aujm — 2a/ auj{auj — m) . 

By returning to the r variable we arrive at the following expressions for the poly- 
nomial solutions of the Heun equation for s = 1: 



(+ii7m(a;;r))i,2 = ^ (-r±-^Jau{auj-m) \ . (3.8) 

r+ — r_ \ UJ J 



Similarly, for s = — 1 we have: 



{^iH.^{uj] r))i,2 = ( r ± -^/auj{auj - m) . (3.9) 

r+ — r_ \ (jj / 

Putting together (13. ip with the specific values of ^, rj, and ( and the polynomial 
solutions of Heun equations, the solutions to the TRE for s = 1 and for s = — 1 can be 
written (modulo normalizing constants) respectively as: 
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{+iRm{uj;r))i^2 = — 7 — (+ii^„(w;r))i,2, (3.10) 

A \ r — r_ 



and 



_,R^{u;■,r))^,2 = e'^''*['-^] (_iif™(u;; r))i,2 , (3.11) 



where 

r^, = r + a+ In |r — r+l — a_ In |r — r_| 

is the "tortoise" coordinate and a± are defined in (13. 5 p . Note that these exact solutions 
for |s| = 1 are presented here for the first time in exphcit form. A similar form of 
polynomial solutions to TRE in the case |s| = 2 can be found in 



Using the orthogonality relations of Heun polynomials [28] it can be shown that in 
terms of the intermediate variables gZ = ±z we have: 



/ A^{sRm{(^; sz))j{sRm{uj; sZ))ld{sZ) = 0, j^l j, / = 1,2, 

J — oo 



(3.12) 



where A is the standard factor from the Kerr metric, defined in (12. 2p . The behavior 
of these solutions at infinity and at the event horizon can be readily determined form 
(I3.10p and (13. lip . First, both solutions in (I3.10p with s = 1 have the behavior: 



r '^e *'^'"* r — > oo (r* — i- oo) 

.iRm{uj;r))i^2 ~ { , (3.13) 

^-Sg-rror. r — ^ r+ (r* — )■ —oo) 

while both solutions in (13. lip with s = — 1 behave like 

j,-(2s+i)gjwr. r OO (r* — )■ oo) 

.ii?^(a;;r))i,2~ <( , (3.14) 

e^'^^* r ^ r+ (r* — > — oo) 

where in the expressions for the behavior at r+ we have set w = u — mQ^. Thus 
we found exact solutions to TRE, the nature of which depends on the relative sign 
between u and w, in agreement with the general analysis in |l]-[6]. When u and tu 
have the same sign then the solutions we found describe one-way waves traveling from 
r+ to infinity or in the opposite direction. The solutions with a; and zj with opposite 
signs describe either waves traveling towards r+ and towards infinity or leaving from 
r+ and coming from infinity. 



9 



3.4 Polynomial solution for spin |s| = 1/2 

This case follows along the same lines like the electromagnetic perturbations but is 
simpler since a = 0. This leads to c^{^ = in (110. Sp which corresponds to a constant 
solution to the Heun equation. This implies that a = and thus determines the 
dependence of the separation constant E from cu. The result is [39l SO] 

:iErr,(aoj) = —a?' Up' + laujm . 

The solutions to TRE for s = ±^ are presented here for the first time: 



.R^{u- r) = i iH^uj) (3.15) 



2 



ii?^(a;;r) = e^-'^*( ^— ^) _iRmioj), (3.16) 



where \Hjn{u) and _iHm{uj) are constants in r. Note that these solutions satisfy the 
same orthogonality relations as those in f l3.12p and have the same behavior at infinity 
and at r+ as those given in ( I3.13P and ( I3.14p . 

4 Solutions to TAE 

4.1 Transformation of TAE into the form of a Heun equation 

In order to transform TAE (12. 5p into the form of a Heun equation we follow the same 
procedure as with the radial one. Following |17]-[5I] we set 

SUoo; u) = {l- + ur^e'^^T^oo; u) , (4.1) 

where u = cos 9 and /xi and /i2 are the indices of the regular singular points at 6* = 
and 9 = 71. We plug (14. ip into (12. 5p and obtain an equation for Tm(w;M). Imposing 
the condition that the equation for Tm(w;M) has the form (110.11) leads to a system of 
quadratic algebraic equations for the parameters /ii, fj,2, and z/. Their solutions are 
given by 

/ii = ± ^ , /i2 = ± ^ , ly = ±au . (4.2) 
For these values of fii, fj,2, and u the equation for Tm{uj; u) is given by 

(PT^{uj] u) 2/ii + 1 ^ 2/i2 + 1 ^ ^ \ dTm{uj]u) ^ 



du^ \u — 1 u + 1 ] du 

aujs^ 



-^-—-{2u U +^2 + 1 + 

1 M+1 V 



u 



(m-i)(m + i)' y ^ p 



10 



In order to complete the transformation we introduce new independent variables (the 
reason for the specific form of this transformation will become clear latter) u i— )■ for 
s = +1/2, 1, 2 and m i— )■ _x for s = —1/2, —1, —2, where ±x = (ld=M)/2. It is important 
to notice that these two independent variables are related by the transformation 6 i— 
71 — 6. Thus we arrive at the following specific non-symmetric canonical form of a 
confluent Heun equation: 

with the identification ±p = ±u, +7 = _(5 = 2fi2 + 1, -7 = +5 = 2/ii + 1, and for both 
signs we have 



aus 

±a = /ii + /i2 + 1 H 

V 



±cr = -2z/(±a) + + a^w^ - {iix + /Ua)^ - (/ii + /ia) - 2z/(/ii - /is) . (4.5) 

The solutions of equation fl4.4p can be expressed with the use of the Frobenius solutions 
a,'y,S,(j; ±x) and all other solutions obtained from it by interchanging the 
finite singular points or performing appropriate 5-homotopic transformations (See the 
Apendix). Thus we obtained the needed solutions to the Heun equation obtained from 
TAE0. In the general case these are given by two sets of infinite series about each one 
of the regular singularities. For more detailed description of all local solutions to the 
angular Teukolsky equation see [391 HO] • With these we can build the solutions to TAE. 

For each choice of a triple /xi, fj,2, and z/ we get a different solution to TAE. There 
are two special cases though: First, we can choose fj,i, fj,2, and z/ in such a way so to 
obtain solutions that are regular at both ^ = and 6 = tt. This chice leads to a well 
studied Sturm-Liouville eigenvalue problem |17]-[5T] for the spin- weighted spheroidal 
wave functions [52j and the separation constant Eim = Eim{auj)^ which for 

fixed s, m, and au are labeled by an additional integer /. The eigenfunctions Simioj] 6) 
are complete and orthogonal on < ^ < vr for each set s, m, and acu. In the case 
s = 0, Sim{oj; 6) are the spheroidal wave functions. When acu = 0, the eigenfunctions 
are the spin-weighted spherical harmonics sY^™" = sS'™(^)e*™''^. In the general case, as it 
is shown for example in [19] the function Tm(w; z) can be expanded as a series of Jacobi 
polynomials in the case of an integer spin, and to their spin-weighted generalizations 
for half- integer spins [U [9l [10], and one can obtain the separation constant A/m as a 
power series in au. 

The second special case corresponds to polynomial solutions to TAE [39| HQ] . Let 
us choose jii = — (m + s)/2,fi2 = {m — s)/2,u = —au. This choice leads in fl4.5p to 
a = 1 — 2s , which for s = 1/2, s = 1 and s = 2 is zero or a negative integer number. 
This means that the condition (110. 9p . which is necessary for having a polynomial 

^The Tome asymptotic solution at infinity Hc^'^\p,a,"/,5,a;x) is not of pliysical significance in 
the case of angular equation (14.31) . because here we are interested only in solutions on the interval 

we (-1,1). 
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solution (a constant in the case a = 0) of the confluent Heun equation is met for these 
values of the spin weight s. The case of zero spin is again excluded by this condition. 
If instead we have s = — 1/2, s = — 1 or s = —2, we must choose /xi = (m + s)/2, fi2 = 

— (m — s)/2, V = au, and this will give us a = 2s + 1, which in this case again is zero 
or a negative integer number. We have to emphasize that this result comes at a price 

- depending on the specific values of m and s either one or the other of the pre-factors 
in (14.11) diverge at the corresponding singular point. This means that for the cases 
of neutrino, electromagnetic and gravitational perturbations we eventually may write 
down the solutions to TAE into the form of diverging at some of the singularities pre- 
factor multiplied by a polynomial expression, simultaneously regular at both regular 
singular points. Again, the neutrino case is relatively simple and can be easily deduced 
from the electromagnetic one, which we will present in detail. 



4.2 Polynomial solutions for perturbations with spin |s| = 1 

Since (110. 9p is satisfied we can obtain a polynomial solution of (14. 4p by imposing as an 
additional (already sufficiency) condition (110. lOp . In our case (110. lOp translates into 
an CO dependence of the separation constant E^, which enters in the parameter a from 

In this case a = —1, so we will have c^'' = in (110. lOp . Because of the specific choice 
we made when introducing ±x in this case we have for both positive and negative spin 
weights the same values of the parameters in the Heun equation: p = —au, a = —1, 7 = 
m,6 = —m, a = Em + a'^u^ + 2au — 2aum. This means that we obtain the same Heun 
equation for both s = 1 and s = — 1 with Frobenius solutions Hc^"-\p, a,'j, S, a; ±x). 
The only difference between the two cases is in the definition of the independent variable 
±x = (1 ± cos^)/2. Thus in this section we will look for a solution only for the case 
s = 1 and will obtain the solution for s = — 1 by performing in the final results the 
transformation 6 ^ tt — 9. In terms of the coefficients from the three term relation 
(110. 6p the sufficiency condition for m 7^ has the form g^^g'f^ = h^f* flf^, where 
g^"* = —a, gl""^ = Aau — a, f^^^ = — m, h^^^ = Aauj. Since this condition leads to a 
quadratic equation for a we obtain pairs of solutions. For both s = +1 and s = — 1 we 
get the same expressions ai = 2au) + 2^yau{au — m) and a2 = 2au — 2^J auj{auj — m). 
Returning back to the original variables of the angular equation we can write down the 
following expressions |37l [3H], for the solutions for s = +1 (for m 7^ 0): 

{+iSm{<^]9))i^2 = — ^-^(cot-) (+iT^(w;6'))i,2, (4.6) 



sin^ V 2^ 

and we have restored the pre-subscripts denoting the two different spin weights. The 
polynomial parts of the solutions for m 7^ are: 



2auj + 2^ au{auj - m) ^6 

_iT„(w;6'))i = 1 ^ cos -, 

m 2 



2au — 2\J au}{ 



au 



m] 



■ cos 



m 
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In the case m = it is easy to solve the equation directly and arrive at the following 
overall solutions: 

UiSoiuj; 6)), = e'"-™^^ tan ^, UiSo{uj; e))2 = e"'^-™^^ cot ^ . 

The behavior of the solutions to TAE we found at the two singularities ^ = and 6 = tt 
can be easily obtained from the expressions above. Each one is divergent either at the 
one or at the other singularity. The exact expressions for m 7^ are: 



^-("^+1) at ^ for m > 1 

^+i^^(a;;^))i,2~ <( . (4.7) 



TT 



ate ^TT for m < 1 



For m = we have {^iSo{u; 9))i ~ 6* ^ at 9 ^ and {^iSo{co; 9))2 ^ {it — 9) ^ at 
9 ^ 7[. If we perform the transformation 9 n — 9 we obtain the solutions for s = —1: 

aw cose / (qr ^ 

(_i5^(a;;^^))i,2 = — ^ cot^— ^ (-iT^(c^; ^))i,2 , (4.8) 

sm V \ z y 

where for m 7^ 0: 



2auj + 2^auj{au-m) . ^9 
(_ir^(w;6'))i = 1 ^- sm - 



2auj - 2^auj{auj-m) . ^ 9 

^ iT^{u;9})2 = I sm - 

m 

and for m = 0: 



USoiou; 9))^ = e'^-=°^^ cot -, USo{co- 9))^ = e^^'^^'tan - . 

It can be shown that because of the orthogonality properties of Heun polynomials [28] , 
the following relations hold for the functions (±15*^(0;; x))j separately for s = 1 and for 
s = -1: 



/ {±iSm{u}]±ix))j{±iSm{uj]±ix))id{±ix) = j^l j,l = 1,2. (4.9) 
Jo 

The second condition (110. lOp for a polynomial solution of the radial equation leads 
to the expressions for the separation "constant" Em for both s = +1 and for s = — 1. 
For any m we have [391 SO] 

±iEm{auj)i = —o?uJ^ + 2aum + 2\/ aujiauj — m) 
±\Em{a''^^2 = —a^oj^ + 2aum — 2a/ au{auj — m) . 
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Surprisingly, these expressions are the same both for the TRE and TAE (see Section 
3.2). 

Hence, the second condition (110. lOp does not produce an additional relation be- 
tween Em and u. The main consequence from this result is that we can express the 
separation constant Em as a function of u but the (complex) frequency itself remains 
unconstrained. As a result we obtain a continuous spectrum in u for the solutions of 
Teukolsky's angular and radial equations (12. 5 p and (12. 6p . 

Thus, the surprising phenomenon of simultaneous fulfilment of the polynomial con- 
ditions for angular and radial equations with spin |s| = 1 is related to the existence of 
the continuous spectrum of Teukolsky Master Equation (12. 3 p in the specific boundary 
problem under consideration. 



4.3 Polynomial solutions for perturbations with spin \s\ = 1/2 

Again we have a = and this determines the dependence of the separation constant 
E from u. The result again is [391 130] 



A_lEm(auj) = —a^'uP' + laujm . 

^2 4 

Like in the electromagnetic case both conditions for polynomiality are the same. Thus 
in this case again we obtain a continuous spectrum. The solutions to TAE for s = ±i, 
presented here for the first time, are: 

cot- ^=-, (4.10) 
^ / V sm d 

.Sm{uj;e)=(t^n-] -^==. (4.11) 

These solutions satisfy the same orthogonality relations as in (14. 9p . The behaviour at 
6 = and at ^ = vr can be easily deduced from (I4.10p and (14. lip . 

5 Nonexistence of gravitational one-way waves of 
continuous spectrum 

The derivation of the polynomial solutions for perturbation with spin |s| = 2 can be 
found in |l5l |l6]. Our analysis in the gravitational case follows the same lines as in 
the electromagnetic one and reproduces the results of these articles. As it was already 
discussed, with appropriate choice of the powers in (13. ip and (14. ip we achieve a = —3 
for both TRE and TAE. When we impose the second polynomiality condition we obtain 
quartic equations for the values of the separation "constant" . The important difference 
is that unlike in the electromagnetic and the neutrino cases, the Heun polynomial 
conditions are different for the angular and for the radial equations. The proof based 
on Taylor series expansions of the roots of quartic equations can be found in [39| HO] . 
This result is consistent with the observation by Teukolsky and Press [6], developed 
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further by Chandrasekar [12j, that the difference between the Starobinsky's constant for 
the angular and for the radial equations in the gravitational case is equal to (12Ma;)^. 
Thus in the gravitational case we have two independent conditions, relating and u 
which leads to a discrete spectrum of u labeled by m and some additional indexes for 
the different solutions of the quartic equations. 

6 Overall solutions to Teukolsky's Master Equation 
for spin |5| = 1 

6.1 Polynomial in both r and cos^, diverging at ^ = and 9 = tt 
solutions 

At this point we are ready to return to the original physical fields. For s = 1 and 
s = — 1 we have respectively \E' = <^o and \E' = p~'^f2- With the solutions we found we 
can write: 



sin 9 



cot 



X : exp In 



and 



((V92)m(w;t,r,6',0))i,2 ~ 



sin 9 



r — Tj 



tan 



[r — la cos I 



ima 

exp I In 

r I — r_ 



r — r I 



(_T„(w; 9))i^2{-Hm{(^] r))i,2 



These expressions can be written in a more compact form if we introduce the Kerr 
coordinates with the relations: 



V = t + r^, 



r+ — r_ 



In 



r — r I 



r 1 — r_ 



In 



r — r J 



Thus the final expressions for the solutions to Teukolsky Master Equation will be 



{{(Po)m{(^;t,r,9,(f))) 



^—iLoV^—aui cos i 



1,2 



(r2 - 2Mr + a^) sin 9 



.WrUTUoo;9))i,2UHU^;r))^,2 

(6.1) 



and 



{{^2)m{0J]t,r,9,(f))) 



1,2 



^—iujU^aw cos d 

(r — ia cos sin ^' 



.Wr{.Tm{uJ■9))^^2{-H^{u;■r)),^2^ 
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where we have introduced the following expressions 

= e*(+^) cot (+6/2) , _W = e^(-^) cot (.6/2) (6.2) 

to denote the stereographic projections of a unit sphere parameterized by angles +0, ^6 = 
9 for s = 1 and _0, ^6 = it — 6 for s = — 1 on the complex planes C^w The first 
formula in (16 .2^ describes stereographic projection from the North pole and the second 
one - from the South one. Note that after the transition from real variables {±6, ±0) 
to the complex one one must introduce an additional phase factor exp {—is±(j)) 
in the spin-weighted spheroidal harmonics, due to the back rotation of the basis (See 
the paper by Goldberg et al. in [8].). In the case of spin 1/2 the introduction of such 
a factor exp {T-^±(f>) is equivalent to a transition in what follows from half-integer to 
integer values of the azimuthal number m and a replacement m — )■ ±1/2. 

The two expressions in (16.11) together with their behavior at r — )■ oo and at the 
regular singularities of TRE and TAE provide us with a basis from which we can build 
solutions with specific boundary conditions. The basis describes wave collimated along 
the poles ^ = and 6 = ir. Depending on the values of the parameters u and m, 
when u and w = u — mfi+ have the same sign, the basis describes waves having 
the same direction of propagation at both r+ and at infinity. Otherwise, when u and 
w = u — mQ^ have opposite signs, the waves have opposite directions of propagation 
at r+ and at infinity. One possible application of this basis is to try to explain the 
Central Engine of the Gamma Ray Bursts (GRB), discussed in [37] and [38] . 



6.2 Polynomial in r regular at ^ = and 9 = tt solutions 

We could combine the polynomial solutions to TRE with the spin-weighted spheroidal 
harmonics representing the regular at ^ = and 6 = n solution. The result will be 
a basis of waves which have the same properties in radial direction as those discussed 
above but will not be collimated along the axes of rotation. Possible application of the 
basis in this form is to study the influence of rotating gravitational field for formation 
and evolution of the Supernovae outbursts [391 SO] ■ 



7 Overall solutions to Teukolsky Master Equation 
for \s\ = 1/2, diverging at ^ = and 6 = n 

Combining the results from (I4.10p . (14. lip . (I3.15p . and (I3.16P we can write down the 
expressions for the neutrino components Xo ^-nd Xi as: 

p—iLoVp—auj cos 6 

{XoUu^;t,r,eA) - -^====j==UWr , (7.1) 
yr"' — 2Mr + a^y/sm 6 

—iu)U pOLO COS 9 

(Xi)™(a;; t, r, <P) ~ , ^^^_wr . 

[r — la cos U) v sm ti 
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Using these solutions, which to the best of our knowledge are published for the first 
time, we will show how we can build regular with respect to the 9 solutions. The 
general form of the one-way solutions, based on (17. ip is 



g -iuiV ^—OLJ cost 



Xo(^;t,r,g,0)= — A^{uj){^Wr. (7.2) 

yr'^ — 2Mr + a^ysmU — 

* ^ m=— oo 

The physical model is determined by the amplitudes Amioj). Physically sound solutions 
correspond to amplitudes Amioj), which lead to a finite result after performing the sum- 
mation. Here we consider only a formal example proving the existence of proper choice 
of the amplitudes giving finite results. Let us write down the part of Xo containing the 
potentially singular factor as 



oo n oo 

sin0)-V2 ^ A^{uj){^Wr = J-{UW\ + UW['') J2 AUio)UWr, (7.3) 



2 

m=— oo 

where \+W\ = cot | > 0. At this point we proceed by choosing in an appropriate 
way functions f{ijj,+W) = J2m=-oo ^mi^) i+W)"^ thus defining the amplitudes. For 
example we can choose f{u,+W) = [+W + {-^W)~^ + const) ^ or more generally 

f{uj,^W) = (^P{u,+W) + Q{u, +W~^)) \ where P and Q are arbitrary polynomials 
of degree not less than one. Then the amplitudes Ajn{oj) are the coefficients in the 
Laurent series of the functions f{u, +W) with respect to +W. With this choice it is 
clear that there are no singularities in the limits \+W\ — )■ and \+W\ — ?► oo. It can be 
shown that choosing the polynomials P and Q properly we preserve the collimation of 
the neutrino waves. 

We can go one step further and formally perform the integration in (12. 4p . Thus we 
arrive at the following general solutions to Teukolsky's Master Equation for neutrino 
waves with spin weights s = 1/2 and s = —1/2 respectively: 



Xo{t,r,9,(f)) 



Xi{t,r,9,(j)) 



Fo{V -ia cos 9, +W) 



Fi{U + ia cos 9,. W) 



\/sm9 

where Fq and Fi are arbitrary functions of their respective variables. The fact that 
arbitrary Fq and Fi satisfy Teukolsky's Master Equation can be verified directly. The 
exact forms of Fq and Fi are to be determined by the specific physical situations. 
The above formal examples demonstrate mathematical technics, which make possible 
the application of singular solutions to Teukolsky's Master Equation for description of 
physical reality. 

8 Conclusion 

In the paper we presented an approach for solving Teukolsky Master Equation based on 
the use of the confiuent Heun equation. After separating the variables we showed that 
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both TRE and TAE can be transformed into the non-symmetric canonical form of the 
confluent Heun equation. The transformation depends on a set of parameters which 
when properly chosen lead to polynomial solutions to both Heun equations related to 
TRE and to TAE. The surprising result (which simply means that there should be 
a deeper physical explanation we do not understand yet) is that for neutrino and for 
electromagnetic perturbations we find solutions which have continuous spectrum, but 
for gravitational perturbations this does not happen. The richness of the results we 
found give us the opportunity to construct different types of solutions in accordance 
with the specific boundary problem we want to study. 

There are many possible directions we intend to pursue. One possibility is to see 
if indeed we can explain some of the basic features of the GRB's and of Supernovae 
using the basis we found. Related problem is to investigate further the stability of Kerr 
black holes. 
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10 Appendix: Heun Equation and Heun Functions 

In this appendix we remind the reader some basic information about the confluent 
Heun equation and its solutions, following the notations of the reference |2B] • 

In subsection 9.1 we give a brief description of the non-symmetrical canonical form 
of confluent Heun equation and its solutions. In subsection 9.2 a basic information 
about local Frobenius and Tome solutions around regular and irregular singular points 
is presented. In subsection 9.3 we remind the basic information about polynomial 
solutions in notations of [28]. In subsection 9.4 we describe the correspondence between 
these notations and the notations, used in the basic articles [25], [26] on the modern 
general theory of all kinds of Heun equations and the properties of their solutions. The 
last notations are used in [29], [30]- [IQ]. At present these conventions become more 
popular, because they are used in the computer package Maple, based on the articles 
|25] . [2B] . This package is still the only one for analytical and numerical computer 
calculations with Heun equations and Heun functions. In Maple's Help one can find an 
available and rich collection of relations and properties of Heun functions of all kinds. 

10.1 Non-symmetrical canonical form of confluent Heun equa- 
tion and its solutions 

The general Heun equation is a second order ODE of Fuchsian type with four regular 
singular points. In the present paper we have to solve the confluent Heun equation 
(CHE) for different cases. It is relatively well studied [23]-[29], but there still exist 
essential gaps in the theory. CHE can be obtained from the general Heun equation by 
coalescing of two of the singular points by redefining certain parameters and taking the 
appropriate limits. In this way two regular singular points coalesce into one irregular 
(in general) point. The solutions of the confluent Heun equation are relatively well- 
studied special functions, already included in modern computer package Maple. These 
functions represent non-trivial generalization of known hypergeometric functions, yet 
have richer properties, because confluent Heun equation has one more singular point 
than the hypergeometric. One of the canonical forms of the confluent Heun equation 
is the so called non-symmetric canonical form [28] : 




The only regular Frobenius' type solution to (110. ip about the regular singular point 
2; = is denoted by Hc''"'\p, a, 7, 5, a; z). It is defined for non-integral values of (1 — 7) 
in the domain |z| < 1 by the condition 



In [28] it is called the "angular solution" of the confluent Heun equation. 



Another solution, is the Tome's type asymptotical solution Hc^^'\p, a, 7, 5, a; z). It 
is defined for complex p = |j9|e"'' in the domain \z\ > 1 by the condition: 






Hc^''\p,a,-f,6,a;0) = 1. 



(10.2) 




lim 2"i7c(")(p,a,7,(5,cr;-|z|e-*'^) = 1. 



(10.3) 
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In [2S] it is called the "radial solution" of the confluent Heun equation. 

Different pairs of local solutions can be constructed using the combinations of 
four known independent transformations of the parameters, which preserve the cho- 
sen canonical form of the Heun Equation. For example by interchanging the regular 
singular points Zi = and Z2 = 1: 

z 1 — z, 
one obtains the following new solutions: 

Hc^°-\-p, a, (5, 7, a + Apa] I - z) 

Hc^''^ {-p, a, 5, 7, a + 4pa; 1 - z) . (10.4) 

All possible sets of local solutions to Regge- Wheeler and Teukolsky equations were 
described for the first time in [391 

10.2 Power-series solutions of the confluent Heun equation 

10.2.1 Taylor series expansion about the regular singularity z = 

If we expand the solution Hc^'^^p, a, 7, 5, cr; z) as a power series 

00 

Hc^^\p,a,j,S,a;z) = J2c^:^^' (10-5) 
then we get a three-term recurrence relation for the coefficients c'^^: 

ft^ctl, + gt^ct^ + ht^c^:l, = (10.6) 
c_i = 0, Co = 1 , 

where 

g(^^ = k{k-4p + -f + S-l)-a 

4") = Ap{k + a - 1) . 

The radius of convergence of the series (110. 5p is equal to unity, which is the distance 
to the next regular singular point [28] . 

10.2.2 Laurent series expansion about the singular point at infinity 

Another power series can be constructed at infinity. In general this series is not con- 
vergent but only asymptotic. For the function Hc^^\p, a,'j, S, a; z) we will have the 
expansion 
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Hc^^\p, a, 7, S, a- z) = ^"^ 4'^^^"'- (10.7) 

k=0 

The three-term recurrence relation for the coefficients cj^^ reads 

ft^ctl + gt^ct^ + ht^ctl, = (10.8) 
c_i = 0, Co = 1 , 

with the following expression^ for fj[\ g^^\ and /i^^^ 

gt^ = {a + k){a + k + Ap--f-5 + l)-a 

ft^ = -Ap{k + 1) 

ht^ = -{k + a-l){a + k-^). 

It is easy to show that in general the series expansion fll0.7p diverges 

10.3 Polynomial solutions of the confluent Heun equation 

Let us consider the case in which the parameter a has a fixed negative integral value 

a = -N, N EN. (10.9) 

In this case the coefficient h^^'^^i vanishes for both expansions above. If we impose in 
addition the second condition that 

c!;i.i = (10.10) 

then the recurrence relation breaks down and we obtain a polynomial of A^-th order 
instead of the infinite series. Since the coefficients g^^^ are linear functions of the 
parameter a the equation c^^^^ = is an algebraic equation of (A^ + l)-th order and 
thus it has {N + 1) zeros ctq, Ui, . . . , uat [281 129] . 

In [39ll40] one can find an explicit representation of the coefficient cj^+i = in form 
of a specific determinant Aat+i. This form is most convenient for practical calculations. 

10.4 Correspondence between the notations of present article 
and the notations, used in computer package Maple 

The computer package Maple uses the conventions of the basic articles [25], [26]. In the 
Maple notation HeunC(a, (3, 7, 6, r], z) for the solution (110.21) the parameters a, /3, 7, 5, 77 

^Note that these expressions are somewhat different from those in [28] . 
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are related in the following way with the parameters of the non-symmetrical canonical 
form of confluent Heun equation [28j, used in the present article, too: 



a 



7^-1 

^ Maple = 4pa - Ml + VMa.i. = 2p7 ^ 

In [29| l39t Ho] a modified Maple-like parametrization of the confluent Heun equation 
is used: 

The equation (110. 12p has a uniform shape. This uniform parametrization simplifies 
the explicit expressions for the coefficients in the series (110. 5p and (110. 7p . For the 
parameters /i and u one obtains the following relations with the parametrization, used 
in present article: 

jJi = CF, ^ + u = ipa. 

The first polynomial condition (110. 9p in Maple notations, as well as in the above 
uniform parametrization, reads: 

5 /3 + 7 

- + ^^ + Ar + l = 0. 
a 2 

It yields discrete values S = —a + 7) + + l) of the Maple parameter 6. Hence, 
the name ^Tv-condition 
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